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Abstract. 

We calculate and resum a perturbative expansion of nuclear enhanced power corrections to the 
structure functions measured in deeply inelastic scattering of leptons on a nuclear target. Our results 
for the Bjorken x-, Q 2 - and A-dependence of nuclear shadowing in F^{x, Q 2 ) and the nuclear 
modifications to (x, Q 2 ), obtained in terms of the QCD factorization approach, are consistent with 
the existing data. We predict the dynamical shadowing from final state interactions in v + A reactions 
for sea and valence quarks in the structure functions F^(x,Q 2 ) and xF^(x,Q 2 ), respectively. In 
p +A collisions we calculate the centrality and rapidity dependent nuclear suppression of single and 
double inclusive hadron production at moderate transverse momenta. 



Dynamical high twist shadowing 

Under the approximation of one-photon exchange, the lepton-hadron DIS cross sec- 
tion dOgh/dxdQ 2 °z L^v W^ v (x, Q 2 ), with Bjorken variable x = Q 2 / (2p ■ q) and virtual 
photon's invariant mass q 2 = —Q 2 . The hadronic tensor can be expressed in terms of 
structure functions based on the polarization states of the exchange virtual photon: 
W^ v (x, Q 2 )=e£ v F T (x, Q 2 ) + £^ v F L (x, Q 2 ).In DIS the exchange photon f of virtuality 
Q 2 and energy v = Q 2 / (Ixm^) probes an effective volume of transverse area l/Q 2 and 
longitudinal extent Azjy x xjy/x, where Azn is the nucleon size, xn = 1/ (2romjv) ~ 0. 1 
and ro ~ 1.2 fm. When Bjorken x <^x/y the lep ton-nucleus DIS covers several nucleons 
in longitudinal direction while it is localized in the transverse plane. 



In the lightcone A + = gauge and the Breit frame we identify the natural short and 
long distance separation of the multiple final state interactions from the propagator 
structure of the struck quark, i(y + /2p + )/(xi — x ± ie) (pole term) and i(xp + /Q 2 )y~ 
(contact term) [1]. The two gluon contact exchange is therefore evaluated in a single 
nucleon state. Resumming the A 1 / 3 -enhanced power corrections we find 01: 

F T \x,Q 2 ) « AFr{ ^ 2{Al Q ^ l \ Q 2 ), (1) 



Q 2 



F^(x,Q 2 ) « AFt LT \x,Q 2 ) + -^Ff(x,Q 2 ). (2) 
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FIGURE 1. Left panel: F2(A)/F2(D) calculation of resummed power corrections versus nuclear A and 
Bjorken-jc fl]. Right panel: FiJSn) / ^(C^show evidence for a power-law in 1/Q 2 behavior consistent 
with the all-twist resummed calculation 1 1]. The bottom right panel illustrates the role of higher twist 
contribution to Fi on the example of R — Ol/ot- 



Here, t, 2 represents the characteristic scale of higer twist per nucleon to 0(a s ): 

e = 3 -^^(p\F 2 (X l )\p), (p\F\X i )\p) = lunj-xG(x,Q 2 ). 
8 x^o 2 

The x- and A-dependence of F 2 (A)/F 2 (D), calculated for £ 2 = 0.09-0.12 GeV 2 , is 
given in the left panel of Fig. 1. Comparison to a leading twist shadowing parameteri- 
zation [2] is also shown. The right panel of Fig. 1 indicates the power law nature of the 
nuclear modification to the structure functions. The physical gluon exchange leads to a 
high twist contribution to the longitudinal structure function F^ and enhances the ratio 
R = Ol/ Ot- We emphasize that both leading twist Jj] and high twist shadowing [1] have 
their origin in the final state coherent scattering. This provides a natural explanation of 
the apparent lack of gluon shadowing in the NLO global analysis [4] which is the only 
one directly sensitive to gluons. 
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FIGURE 2. Left panel: power corrections to the structure functions /^(x, Q 2 ) and xF^(x,Q 2 ) |5] for 
two values of xb corresponding to NuTeV measurements 0. Right panel: high twist modification to the 
Gross-Llewellyn-Smith sum rule A GLS \5]. 



Neutrino-nucleus scattering 

Neutrino-nucleus scattering provides the unique opportunity to separately study the 
effect of coherent power corrections for sea and valance quarks [5] through the structure 
functions: 



F$ ) (x B ,(?) » A(2) 



v , . , £ 2 (A 1 / 3 — 1) M 2 , 
£ \V DU \ 2 ^ ( XB+XB^-g L+ XB -JjL^ 

v , , , , £ 2 (AV3_i) Ml 
H-)£\VuD\ 2 $[xB+x B ^-g 2 ; -+x B -^,Q 2 



• (3) 



Here Vdu are the CKM matrix elements. Eq. © identifies the nuclear enhanced high 
twist corrections with dynamical mass m 2 , = <^ 2 (A 1//3 — 1) generated by the final state 
parton scattering through direct comparison to My D . 

The modification to the structure functions F2(x,Q 2 ) and xF^,{x, Q 2 ) for two select 
values of xb are shown in the left panel of Fig. 2. These give a good description of the 
observed power law deviation of the reduced cross sections measured by NuTeV |0, 
0] from the leading twist pQCD at small values of Q 2 . Note the difference in the 
"shadowing" of F 2 and xFt, due to the different steepness of sea and valence quark PDFs 
(in x). The right panel of Fig. 2 demonstrates the improved agreement between data and 
theory for the Gross-Llewellyn-Smith sum rule H: 



Agls = \ (3 - 5 GLS ) = + ^ + m-") • (4) 
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FIGURE 3. Left panel: upper limit on the suppression of the single inclusive particle production 
Rpl (PT\ ) fr° m coherent power corrections versus rapidity and centrality 01 . Data is from BRAHMS 

Right panel: suppression of the double inclusive cross section R p ^(pri ,Pt 2 ) for different rapidity gaps, 
PTdPt 2 ranges and centrality. 



Proton-nucleus collisions 



The p +A analogue of the DIS coherent power corrections is the final state interactions 
of the small x b parton in the \i\ <C \s\, \u\ regime. Here t = q 2 = (x a P a — Pc/zi) 2 and the 
x b rescaling in the lowest order pQCD formalism reads M\- 



ab—>cd 



( x b) 



Fabled (x b l+Q^(AV3_i) \ 



(5) 



In Eq.© F ab ^ cd (x h ) = \M ab ^ cd \ 2 ^)(x b )/x b and C d is a color factor, C q ^ = 1 and 
C g = Ca/Cf = 9/4 for quark (antiquark) and gluon, respectively. 

The left panel of Fig. 3 shows the upper limit on the centrality and rapidity depen- 
dent suppression Rj^ of single inclusive hadron production at RHIC. Data is from 
BRAHMS [9]. Additional nuclear suppression arises form the energy loss in cold nu- 
clei iflpll . The right panel shows the suppression of away side dihadron correlations R^ 
versus transverse momentum, rapidity and centrality on the example of the area of the 
correlation function C(A0) = dN hl,fl2 / dA(j) The pronounced pr 2 dependence is consis- 
tent with STAR preliminary data ifTTIl . 
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